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Abstract 

We study the transport properties of passive inertial particles in a 2 — d incompressible flows. 
Here, the particle dynamics is represented by the A — d dissipative embedding map of the 2 — d 
area-preserving standard map which models the incompressible flow. The system is a model for 
impurity dynamics in a fluid and is characterized by two parameters, the inertia parameter a, and 
the dissipation parameter 7. The aerosol regime, where the particles are denser than the fluid, and 
the bubble regime where they are less dense than the fluid, correspond to the parameter regimes 
a > 1, and a < 1 respectively. Earlier studies of this system show a rich phase diagram with 
dynamical regimes corresponding to periodic orbits, chaotic structures and mixed regimes. We 
obtain the statistical characterizers of transport for this system in these dynamical regimes. These 
are the recurrence time statistics, the diffusion constant, and the distribution of jump lengths. 
The recurrence time distribution shows a power law tail in the dynamical regimes where there is 
preferential concentration of particles in sticky regions of the phase space, and an exponential decay 
in mixing regimes. The diffusion constant shows behaviour of three types - normal, subdiffusive and 
super diffusive, depending on the parameter regimes. Phase diagrams of the system are constructed 
to differentiate different types of diffusion behaviour, as well as the behaviour of the absolute drift. 
We correlate the dynamical regimes seen for the system at different parameter values with the 
transport properties observed at these regimes, and in the behaviour of the transients. This 
system also shows the existence of a crisis and unstable dimension variability at certain parameter 
values. The signature of the unstable dimension variability is seen in the statistical characterizers 
of transport. We discuss the implications of our results for realistic systems. 
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I. INTRODUCTION 



The study of transport properties of impurities in fluid flows is a problem of important 
practical interest. These properties can have serious implications for pollution in the atmo- 
sphere, plankton populations in the ocean, and diverse engineering applications. There is 
evidence that transport processes occurring in nature, e.g. the atmosphere 111, are quasi- 
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two dimensional [2|, |3|]. The chaotic advection |5|] of impurities which are modelled by 
of inertial particles of finite size in flows of different types has been studied both in the 
case of passive particles , and in the case of active particles that can react with 

the surroundings j^, [l^. Impurity dynamics in such two dimensional fluid flows can be 
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effectively modelled by the bail-out embedding maps of 2 — area-preserving maps 
We attempt to understand the dynamical and statistical properties of impurity transport in 
2 — d incompressible flows by studying such embedding maps. 

The Lagrangian dynamics of small spherical tracers in two dimensional incompressible 
fluid flows is described by the Maxey-Riley equations. These are further simplified under 
various approximations to give a set of minimal equations called the embedding equations 
where the fluid flow dynamics is embedded in a larger set of equations which include the 
the differences between the particle and fluid velocities [ll|, [l^. Although the Lagrangian 
dynamics of the underlying fluid flow is incompressible, the particle motion is compressible 



14( 1 . and has regions of contraction and expansion. The density grows in the former giving 



rise to clusters and falls in the latter giving rise to voids. The properties of the base flow 
have important consequences for the transport and mixing of particles. Map analogs of 
the embedding equations have also been constructed for cases where the fluid dynamics is 
modelled by area-preserving maps which essentially retain the qualitative features of the 
flow [15]. The embedded dynamics in both cases is dissipative in nature. Here, we study 
the embedded standard map. 

We study several statistical characterizers of transport properties in different dynamical 
regimes of the system. These are the recurrence time statistics, diffusion and drift quan- 
tiflers, and the distribution of jump lengths. The flrst of these, viz. the recurrence time 
statistics, show signatures of the dynamical regime of the system. The cumulative proba- 
bility distribution of the recurrence times (sometimes called the complementary cumulative 
distribution [16]), in the regime where chaotic structures are seen shows the power law tail 
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characteristic of dynamics in an inhomogeneous phase space with sticky regions. In contrast, 
in the regimes where full mixing is seen the RTD shows exponential decays characteristic 
of rapid mixing. Earlier studies of the connection between the recurrence time statistics 
and the transport of particles had established that normal transport occurs for cases where 
the recurrence times distribution shows exponential behaviour, while anomalous transport 



occurs if there is a power law recurrence time distribution [2J, |25|] . In the embedding map 
studied here, we observe that there is a crossover from exponential behaviour to power law 
behaviour in the recurrence time distribution in different dynamical regimes implying that 
there are some dynamical regimes which are characterized by anomalous transport. The 
diffusion properties corresponding to these cases are further explored. We also comment on 
the differences between the behaviour of aerosols and that of bubbles. 

The diffusion studies show that the system can show a variety of regimes corresponding to 
generalized diffusive behaviour. Three types of diffusion, viz. normal or brownian diffusion, 
superdiffusion, and subdiffusive behaviour are identified. A phase diagram is constructed to 
classify and demarcate the regions of diffusion in the a — 7 parameter space. Apart from the 
usually observed normal and superdiffiusive cases, the embedding map has two subclasses 
of subdiffusion - one associated with trapping with stationary states and the other with 
trapping with non- stationary states. The regions where ballistic diffusion is observed in the 
phase space are identified in the a — 7 phase diagram. Similarly, the areas where the average 
drift is zero are also likewise identified in the phase diagram. The diffusive behaviour of the 
system is compared with the dynamical regimes. Our inferences can have implication for 
the transport of pollutants in the atmosphere or other real application contexts. 

The signature of the different dynamical regimes can also be seen in the transient be- 
haviour of the system. The jump length distributions, where the jump length is defined as 
the Euclidean distance between successive iterates, have envelopes which can be fitted by 
a heavy tailed Levy distribution in the chaotic structure regime which are characterized by 
sticky regions in the inhomogeneous phase space, and by a Gaussian envelope in the mixing 
regimes. The cumulated jump length distributions here show power-law behaviour for the 
sticky regimes, and cumulated Gaussians in the mixing regimes. 

In the aerosol regime of the embedding map, the plot of the largest Lyapunov exponent 
shows signatures of unstable dimension variability (UDV) at certain parameter values. The 
recurrence time statistics shows the signature of the UDV, as does the diffusion constant. 
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We discuss the consequences of this UDV for the transport properties of the system. We 
summarize the imphcations of our results. 



II. DYNAMICAL REGIMES AND PHASE DIAGRAM 



The dynamics of inertial particles in a flow has been shown to be modeled by the bailout 



embedding equation 



dv du 

Here, the velocity of the particle is v = d'^/dt and the velocity field of the fluid flow is 
u(x, for a 2-d fluid. The inertial parameter, a is related to the particle and fluid 
densities, pp and pj respectively, by the equation a = 3pf/{pf + 2pp). Thus, the a < 1 
regime corresponds to the aerosols, and the a > 1 regime corresponds to the bubbles. The 
dissipation parameter 7 is defined by 7 = 2a/3St where St is the Stokes number, and 
provides a measure of the contraction or expansion in the phase space of the particle. The 
^article dynamics is dissipative in nature. A map analog of Eq. [1] has been constructed 
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13|. This has the form 



X„+i = M(X„) + 5n 

Sn+i = e-^[ax„+i-M(xJ]. (2) 

Here, the area preserving map M(x) represents the incompressible fluid acting as the base 
flow. The vector x represents the position of the particle, and the vector 6 defines the 



detachment of the particle from the fluid velocities [ll|. In our study, the dynamics of the 
base fluid flow has been represented by the 2 — d standard map, as it is a prototypical area 
preserving system j^^], and is widely used as a test bed in a variety of transport problems 



281]. The evolution equations of the standard map are. 



K , 

Xn+\ = a;„ + — sm(27r?/„), 

Vn+l = yn + Xn+1 (3) 

with periodic boundary conditions such that x,y ^ [—1, 1]. We fix the parameter K = 2 
in our study. The phase space of the standard map contains both regular and chaotic 



regions for this value. Substituting the form of the standard map as M in the embedding 
map equation above, we get a 4-dimensional embedding map, 



K 

Xn+i = Xn + — sin(27r?/„) + 



Vn+l 



K 

Xn + yn + ^ sin(27r?/„) + 51 
Zn 
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n+1 



iVn+l - 



(4) 



This map is invertible and dissipative. We study the behaviour of the map in the pa- 
rameter regimes < a < 3, and < 7 < 1. Thus, our study encompasses both the aerosol 
regime i.e. a < 1 and the bubble regime a > 1. The dynamical regimes of the system have 
been catalogued in the phase diagram of Ref. 12| for the same range of parameters. A 
variety of dynamical behaviours can be seen at different values of a and 7. These include 
periodic behaviour, chaotic structures (Fig. [3](a)) and fully mixing regimes Fig. [5]^a). A 
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detailed phase diagram of the system can be found in Fig{T] ( See Ref. 
description of the method of construction of the phase diagram) |40l ]. 

Since the map represents the Lagrangian dynamics of particles, the chaotic structures 
seen due to the stickiness of the phase space correspond to the preferential concentration of 
particles. On the other hand, a uniform distribution of particles can be seen in the well mixed 
regimes of phase space. Thus, the phase diagram of the system has important implications 
for the transport properties of the system. In the next section, we examine these transport 
properties via statistical characterizers. We start with the recurrence time statistics. 



III. RECURRENCE TIME STATISTICS 



The statistics of recurrence times are of fundamental importance in the study of chaotic 



systems 
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These involve the study of recurrences of a given dynamical state of the 
system in finite time. The first recurrence time of a trajectory of a system can be defined 
to be the time ri taken for a trajectory, which starts from a small subset ^ of the phase 
space of the system, to return to the same subset ^ in the limit where the volume of the 
subset ^ — > 0. The trajectory may return to the subset at subsequent times T2, T3, ...Tn with 
Tn being called the n^^ recurrence time. A set of the time intervals Tj can be obtained in 




FIG. 1: (Color online) The phase diagram showing the three dynamical regimes in the embedding 
map in the a-7 parameter space. The regions marked by red ('M') show mixing and the regions 
marked by blue ('P') have periodic orbits. Chaotic structures are seen in the regions marked by 
green and yellow ('C'), in the aerosol and bubble regions respectively. There are some regions 
marked by light blue (labeled by '*') which have chaotic structure over a mixing background in 
phase space 

the long time limit to give the recurrence time distribution of that trajectory. The average 
recurrence time of the subset ^ is calculated by averaging over the recurrence times of the 
trajectories starting in the subset ^, and the average recurrence time of the entire phase 



space can be obtained 
phase space. See Refs. 

Extensive studies o 
Hamiltonian systems 
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averaging over the recurrence times of all the partitions in the 



26| for more rigorous definitions. 
' the mean recurrence times have been carried out in the case of chaotic 
32| . Generally, Hamiltonian systems are not fully hyperbolic, and the 
phase space of such Hamiltonian as well as that of area preserving systems contains regular 
islands and chaotic mixing regions in the phase space. A trajectory that originates in the 
chaotic region can stick to the neighborhood of the regular islands intermittently, and as a 
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FIG. 2: The schematic illustration of the recurrence phenomena in an invariant set T. A trajectory 
starting in a small partition ^ is revisiting the partition in finite time. Here only the first and the 
second recurrences are shown. 



result, influence long time properties like the recurrence times giving rise to power law decays 



in the recurrence time distribution 



33| . The cumulative recurrence time distribution shows 



exponential decay for uniformly hyperbolic system, where mixing is strong. In contrast, 
power law decay is observed for the recurrence time distributions of systems that have 
inhomogeneous phase space [34 1. 

In the case of area preserving maps, such as the standard map, the Poincare recurrence 
time statistics shows two limiting cases. Firstly, for the chaotic strong coupling limit (high 
values of the nonlinearity parameter K), where mixing is prominent in the phase space, the 
recurrence time distributions show exponential decay as functions of time. Secondly, for the 
near integrable weak coupling limit (small values of K, K < Kcruicai = 0.97163540631), the 
distributions show inverse power law behaviour [35|. At K = KcnUcai the distribution shows 
a power law decay [32|]. 

The base flow of our embedding map is the area-preserving standard map, with the 
recurrence time statistics described above. However, the dynamics of the inertial particles, 
as described by the embedding map, is dissipative in nature. It is therefore interesting to 
study the recurrence time statistics of the embedding map and see how the statistics of 
the inertial particles differs from the statistics of the base flow, and also see the effects of 
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FIG. 3: (a) The phase space plot with chaotic structures (a = 0.4 and 7 = 0.4). The darker regions 
are visited more often and the recurrence times in these areas are longer, (b) The 3-dimensional 
invariant density plot of the attractor in the x — y plane shows larger values of the invariant density 
p where the region is sticky, (c) The histogram is plotted to show the distribution of the normalized 
invariant density p. 

different dynamical regimes. We do this in the current section. 

We study the recurrence time statistics in the 2-dimensional phase space spanned by the 
configuration space coordinates x — y and compare this with the recurrence time statistics 
in the 4-dimensional phase space {x, y, Sx, Sy). For the 2 — d statistics, the two dimensional 
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FIG. 4: (a) The complementary cumulative distribution of recurrence times at this parameter 
value shows a power law decay in the long time limit, with the exponent (^2d = Cm ~ ^ — ~2.9 for 
the 2-d recurrence and ^4,^ = — 1 = —2.6 for the 4-recurrence. (b) The log-linear plot in the 
short time limit shows an an exponential decay, with the slope = —0.00134 for the 2-d case and 
slope = —0.0005 for the 4-d case. 

configuration space x, y with x G [—1,1] and y G [— 1, 1] is divided into a grid of 50 x 50 boxes 
and 200 uniform random initial conditions are evolved for 1, 000, 000 time steps, with 5, 000 
iterations as transients. Each trajectory is associated with the box it visits immediately 
after the transients. The average recurrence time associated with any box is the recurrence 



time averaged over the initial conditions that are associated with that box j4l| . We compare 
the general features of the recurrence time distributions for the two aperiodic dynamical 
regimes, viz., the chaotic structure regime and the mixing regime (see Fig. [1]) . 

First, let us consider the chaotic structure regime with the parameters a = OA and 
7 = 0.4. We first discuss the 2 — d recurrences. It is clear from the x — y phase space plot 
of Fig. ISt^a) that the phase space in the chaotic structure regime is inhomogeneous. If the 
region available to all the trajectories for motion is called F, there are some special regions 
5*, contained in F where the trajectory visits more often than others. Hence, in the phase 
space plots, the region S is darker than the background F (see Fig. [3]^a)). These darker 
regions are identified as the sticky regions, and have an important effect on the recurrence 
time statistics. The invariant density p supported by the x — y space is plotted in Fig. 



[3]^b). To obtain the invariant density, random initial conditions were uniformly spread in 
the phase space covered by a grid, and the total number of times that trajectories visit a 
particular box in its itinerary was counted 2l|. It is clear that the value of the invariant 
density in the sticky regions is much higher than that in the background. The histogram 
of the normalized invariant density is plotted in Fig. ^c) and clearly shows that the major 
part of the phase space (the background) supports very small values of the invariant density, 
and the large values of the invariant density are concentrated on a small subset of the phase 
space, the sfcky .eg.,., or the regions where the ehaotie structure .s seen Q. Such sticky 
regions have also been found in area preserving maps and experiments in fluids [6|, |22|. 

The cumulative recurrence time statistics in the chaotic structure case show a power law 
decay for large times (see Fig. 111(a)), 



P(r) ~ T<'. (5) 

Here, the exponent of the cumulative probability distribution and the exponent ( of 
the probability distribution are related by ( =C — 1, as both the distributions show power 
law behaviour at large times. In contrast, at short times, the cumulative RTD shows an 
exponential decay (see Fig. 111(b)), 



P(r) ~ e^. (6) 

We compare the 2 — d recurrences seen here with the A — d recurrences in the full x — y, 
6x — 5y space. The 4-dimensional recurrence time distribution corresponding to the chaotic 
structure case is plotted along with the 2-dimensional recurrence time distribution in Figs. 
m^a)-(b). Since the 2 — d space is a projection of the A — d space on the x—y plane, recurrences 
to the A — d space are expected to be longer than the recurrences to the 2 — d space, hence 
the distribution in 4 — has shifted to higher r values. The character of the graph is 
unaltered, namely it shows a power law decay though the exponent is = —2.6, while the 
2-d exponent is = —2.9. In the short time limit, where the distribution is exponential 
for both cases, and the constant h takes the value = —0.00134 and hid = —0.0005 for 
the 2-dimensional and the 4-dimensional recurrences respectively. We note that the values 
of ( and b are not universal and change with the parameter values of a and 7. However, 
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the qualitative behaviour of the recurrence time distribution remains the same at other 
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parameter values in the chaotic structure regime 

In the case of the area preserving maps, the power law scaling indicates that the phase 
space is inhomogeneous. In the case of the dissipative map discussed here, a similar effect is 
seen due to the existence of the chaotic structure or sticky regions in the phase space. The 
examination of the residence times of trajectories starting from arbitrary initial conditions, 
shows that the trajectories spend 80% of their time in the sticky regions, regardless of the 
initial conditions. Trajectories which start in the non-sticky regions of the phase space wan- 
der into the sticky regions and spend a long time there before emerging outside. This leads 
to long recurrence times. Other trajectories have short recurrence times. The exponential 
decay seen in the early part of the RTD corresponds to the short recurrences. On the other 
hand, the power law scaling regime is due to the recurrences to the non-sticky regions. This 
accounts for the crossover from exponential to power law behaviour seen in the recurrence 
time distribution. Thus, the chaotic structures play the same role in the phase space for 
the embedding map as the islands seen in the area-preserving maps. However, there is no 
hierarchy in the chaotic structures, unlike that in the case of the hierarchical islands of the 
area-preserving maps. 

The behaviour of the recurrence time distribution in the mixing regimes of the phase 
diagram is quite different. Here, unlike the chaotic structure case, the mixing is very strong 
and there are no preferential regions in the phase space which the particles visit more 
frequently, or stick to for longer times. The phase space plot of a typical mixing regime 
can be seen in Fig. [5](a). The recurrence time distributions show exponential decay both 
in the 2-d and the 4-d CcLSG, cLS IS expected for the strongly mixing case. In Figs. [5] (c) 
and (d) the log-linear plots of the RTD shows the exponential decays with the slopes for 
the 2-dimensional and 4-dimensional cases to be —0.000175 and —0.0000011 respectively. 
The histograms of the recurrence times for the 2-dimensional and the 4-dimensional cases 
are shown in Fig. [5]^b). It is clear that the 4-dimensional recurrence times are longer than 
the 2-dimensional recurrence times. Also it is evident that the frequency of occurrence of 
the 4-d recurrences is much smaller than the frequency of occurrence of the 2-d recurrences 
(compare the values of /(r) in both the ordinates of Fig. [5](b)). 

In the chaotic structure case we saw a crossover from the exponential decay behaviour at 
short times to power law behaviour at asymptotic times. This crossover behaviour seen in 
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FIG. 5: (a) The phase space plot in a mixing regime (a = 2.8 and 7 = 0.4). We can see there is 
no preferential concentration of particles in any region of the phase space, (b) The histograms of 
the distribution of recurrence times in 4-dimensional case (top) and 2-dimensional case (bottom). 
The log-linear plots of the complementary cumulative distribution of recurrence times in (c) 2- 
dimensions with slope —0.000175, and (d) 4-dimensions with slope —0.0000011. 
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the recurrence statistics signals a transition in the transport properties of the system 
We examine the transport properties using detailed diffusion studies, and correlate them 
with the recurrence times and the dynamical regimes. 
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IV. DIFFUSION 



The transport of passive inertial particles in flows can be described statistically by ex- 
amining the dispersion as a function of parameters a and 7. We consider an ensemble of N 
initial particles distributed uniformly randomly in the phase space and evolve them in time. 
As the initial conditions evolve in time, the particle cloud drifts in the two dimensional con- 
flgurational space from the initial position as well as the individual particles disperse from 
the moving cloud. The dispersion of particles is given by the variance of the displacement 
of particles a^, 



Here x(t) is the position of a particle and < x(t) > is the average position of all the 
particles at time t, both in the conflguration space. The diffusion coefficient D and the 
exponent 77 quantify the type of diffusion. The angular brackets indicate the average over the 
ensemble. The conflguration space in two dimensions considered here is the cover space, i.e., 
without the periodic boundary condition that was used in the previous section. Generally, 
when the exponent in Eqn. [7] takes values r/ > 1, the process is called superdiffusion and 
the trajectories of the particles have long displacements. The transport is characterized by 
normal diffusion if the variance grows linearly with time, i.e. if 77 = 1. Subdiffusive transport 
occurs when the exponent 77 takes values less than 1. The embedding map shows all the 
three main classes of diffusion processes as illustrated in Fig. [61 Subdiffusive behaviour in 
the embedding map can be further classifled into two subclasses, viz., one associated with 
the trapping regions with non-stationary states and the other with the trapping regions with 
stationary states. We analyze the behaviour of these classes in detail below. 

The phase diagram of Fig. [1] separates out the distinct dynamical regimes in the a — 7 
space. Here, we flnd it useful to obtain a phase diagram which classifles the main diffusion 
regimes as a function of the parameters a — 7. The value of the exponent r/ can be used for 
such a classiflcation. For this, the log-log plot of the variance as a function of time is fltted 
to a straight line after the initial transients. We used a linear square flt for flnding the value 
of ?7 for each data point in the a — 7 parameter space. 

Theoretically, the normal diffusion regions in the phase diagram can be distinguished 
from the anomalous diffusion regions if the value of the exponent r/ > 1 and < 1. When 
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FIG. 6: The plot showing the variance of the particle cloud as it evolves in time, (a) Superdiffusive 
behaviour is observed for parameters a = 1.2 and 7 = 0.2 with exponent rj = 2.0135. (b) Normal 
diffusion is seen for the parameters a = 2.8 and 7 = 0.1 with the exponent rj = 1.00611 (error 
bars). Trapping regimes show the plateauing of the variance (c) with nonstationary states (q = 0.5 
and 7 = 0.5) and (d) trapping with stationary states ( a = 1.5 and 7 = 0.8). 

dealing with numerical data, the normal diffusion regions cannot be identified by selecting 
the ones with the exponent rj = 1. The distribution of the values of rj in the entire phase 
diagram is shown in Fig. El^a), and shows three peaks, centered around ~ 0, 77 ~ 1 and 
7] ~ 2. The inset shows the distribution centered around 17 ~ 1 that can be fitted with a 
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FIG. 7: The probability distribution of the exponents rj in the phase diagram, showing three clear 
peaks. The inset shows the peak centered around r] = 1 for which the full width at half maximum 
ujjn is in the range 0.835 < r] < 1.165 that will delimit the normal diffusion exponents. The peak 
value at rj = is f{r]) = 6000. (b) The probability distribution of drifts in the phase diagram 
has three classes, viz., high drifts {Idrift > 12), medium drifts (12 > Idrift > 3), and low drifts 
(3 > Idrift) (see inset). 

Gaussian. The standard deviation of the Gaussian is a gauss = 0.14, and the full width at half 
maximum {uim) for this Gaussian, expressed as lu^ = 2(7gauss^/'^log{2) = 2.35A8agauss, turns 
out to be 0.33 (see the inset of the Fig. [7](a)). This Um centered around r] = 1 can therefore 
be taken as a good practical bound on the rj, to demarcate normal and anomalous diffusion 
regimes. The upper and lower bounds on the normal diffusion region are rjiower = 1 — 
and Tjupper = 1 + oJm/2 respectively. The exponent values that demarcate the regions of the 
normal diffusion are therefore rjnormai £ [viower, Vupper]- The cases which have values less than 
the lower bound rjio^^r are identified as subdiffusive and the ones which have values higher 
than the upper bound r/npper are identified as super diffusive. 

We note that in generating the phase diagram for diffusion regimes we have used 3000 
time steps with 100 initial conditions distributed uniformly randomly in the phase space. 
For the phase diagrams in Fig. [H]a resolution of 0.01 in both the directions was used. 

The diffusive phase diagram indicates the following transport properties. The phase 
diagram shown in Fig. [8]^a) has normal diffusion regions in the regions that are labelled 
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FIG. 8: (Color online) (a) The phase diagrams showing the four diffusion regions in the embedding 
map in the Q-7 parameter space, superdiffusive regime is marked by 'S' (blue), normal diffusive 
regimes are marked by 'N' (red), subdiffusion with stationary states are marked by 'T*' (pink) and 
subdiffusion with nonstationary states are marked by 'T' (green), (b) The phase diagram shows 
the drifts of particles in the 0-7 parameter space. The regions marked in 'blue' have large average 
drifts, Idrift > 12 (labelled by 'C'), regions in 'red' have medium average drifts, 12 > Idrift > 3 
(labelled by 'B'), and regions in 'green' have low drifts, Idrift < 3 (labelled by 'A'). This is for 3000 
time steps averaged over 100 random initial conditions uniformly distributed in phase space. 

'N' (colored red), which are predominantly seen in the areas where the mixing is seen in 
the phase space (see Fig. [1]). The superdiffusive regions are labelled 'S' (coloured blue), 
as mentioned above, are identified where rj > 1.165. On comparing these regions with the 
phase diagram of the dynamical regimes, it is clear that the superdiffusive regions coincide 
largely with the periodic regime. 

The subdiffusive region is delimited by the criterion rj < 0.835. It is easy to see that 
these regions mostly coincide with the chaotic structures regimes seen in Fig. [1] Further, 
there are two types of subdiffusion occurring here, viz. those associated with trapping in 
stationary states and trapping in non- stationary states. On the aerosol side of the phase 
diagram (colored green and labelled 'T'), particles behave as though they are trapped in 
an attractor, but their dispersion grows sublinearly with time (Fig. [6](c)). The cover phase 
space for such a typical case is seen in Fig. [TO](a). On the other hand, in the bubble 
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region (colored pink and labelled 'T*'), particle trajectories get trapped forever and become 
stationary. Fig. [TOT b) shows the cover phase space for the case for trapping into stationary 
states. As a result of the motion tending to fixed points, the variance tends to a constant 
value without any fluctuation after the transient time as seen in Fig. E^d). 

There are certain parameter regimes at which the particles do not diffuse beyond a given 
region in the phase space and we saw that they belong to the subdiffusive regions. The 
average drift of the particles can be used to quantify the overall drift of the particle cloud 
from the initial position it started with. The drift is defined as the mean position of the 



cloud of the particles after a time t, i.e., Idrift = l{t) — 1(0), i.e. where Idrift is the difference 
between the average position l[t) of the cloud of the particles in cover space, at time t from 
that at the time t = 0. The distribution of the average drifts of the particle cloud for the 
range of values of a — 7 discussed in the phase diagram, is plotted in Fig. [7](b). It is clear 
from Fig. El^b) that there are three main ranges of drifts. The large drifts are identified by 



the criterion Idrift > 12, the intermediate ranges of drift by 3 < Idrift < 12, while the small 



drifts are identified to have drift rates Idrift < 3 (see inset of Fig. [3(b)). 

A phase diagram for the drifts in the a — 7 space is plotted in Fig. [H](b). Regions of large 
drift are coloured blue, those of intermediate drift are coloured red, and green regions have 
small drift. On comparing the parameter regions 'T*' (pink) in Fig. [8]^a) where particles are 
trapped and become stationary, with the corresponding region of Fig. Mjo), it is clear that 
the particles here have low values of drift. The regions of high drifts in Fig. [8](b) overlap 
with the super diffusion regions 'S' (blue) both in the aerosol and the bubble regime in Fig. 
IHI^a). Therefore, there is good correlation between the phase diagrams constructed from the 
dispersion and drift criteria. 

The phase diagrams obtained in Fig. [8] show no clear contiguous areas defining the three 
diffusion regimes, or the three drift regimes, and these regions are seen to be interspersed in 
some areas of the phase diagrams. This is due to the fact that the values of the slope r] and 
the average drift Idrift belong to continuous distributions which are partitioned into three 
regions. Hence sharp demarcations are not achieved. 

We note that the embedding map also shows a region corresponding to the ballistic 
super diffusion regime 25| which is characterized by a diffusion exponent ~ 2 and has been 
earlier seen in diverse contexts such as the motion of atoms, molecules and clusters on solid 

regime is shown 
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FIG. 9: (a) The phase diagram shows the regions characterized by ballistic diffusion. Here the 
diffusion exponent take values r] = 2. (b) The phase diagram with points marked showing regions 
where the average drift Idrift ~ 0. 

in the phase diagram Fig. [^a). We note that this regime is seen at the edges of the periodic 
tongues. It is also useful to have a phase diagram to locate the regions where the average 
drift is zero, Idrift ~ 0. The Fig. ^Hjo) shows the regions of near zero drift and in reference 
to Fig. Mjo), we can see that it forms a subset of the regimes of the slowest drifts, marked 
in green (labelled 'A'). 

It is important to note that characterization and classification of different diffusion and 
drift regimes are very useful in real application contexts. We expand on this point in the 
conclusion. 

V. THE BEHAVIOUR OF THE INITIAL TRANSIENT 

The phase diagrams of the previous section have been obtained after the systems has 
equilibrated after an initial transient. However, the behaviour of the initial transient can 
itself be quite different in the case of distinct dynamical regimes. We illustrate this for the 
dynamical regimes of the embedding map. 

We examine the probability distribution of the jump lengths, i.e. the Euclidean distance 
between successive iterates, for transient times upto 500 iterates in regimes corresponding to 
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FIG. 10: The phase space plots of the subdiffusion regimes (a) nonstationary states (b) Trapping 
states with stationary points 

periodic, chaotic and mixing behaviour. Here, 1000 initial conditions are spread uniformly 
randomly in the phase space for this analysis. Fig. fTDfa) shows the probability distributions 
of the jump lengths of the trajectories in the regime {a = 0.5,7 = 0.5), where chaotic 
structures are seen in the phase space. The distribution is clearly characterized by a heavy 
tail. Similar heavy tails characterize the jump length distribution of the trajectories seen 
at other points in the chaotic structure regime. In the region of the heavy tail the envelope 
can be fitted by the function. 



flevyi^) 




1 



-exp 



(8) 



2vr (/_ 5^)3/2 --^2(/- 50 
where, the parameters take the values v = 0.42 and 6i = 0.34 and the scale factor is 
K = 11 for the case shown in Fig. [TTT a). On the other hand, the envelope of the distribution 
of jump lengths of trajectories for parameter values in the mixing regime a = 2.8, 7 = 0.1 
shows Gaussian behaviour. In the Fig. [TTT b) the envelope of the distribution has been fitted 



with a Gaussian function, fgaussi}) = )^ exp( 



), with fx = 0.95, a = 0.58. 



Thus, in the initial transient in the chaotic structure regime, the distribution of jump 
lengths peaks at short values, but decays with a power-law tail, characteristic of the sticky 
regions. On the other hand, the jump length distribution in the mixing regime, shows normal 
behaviour, indicating rapid mixing of the initial conditions. 
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(a) / (b) I 

FIG. 11: The envelope of the distribution of flight lengths of trajectories in the periodic and 
chaotic structure regimes is fitted with a heavy tailed distribution, (a) the periodic regime (a = 
0.5,7 = 0.5). (b) In the mixing regime the envelope is fitted by a Gaussian (a = 2.8,7 = 0.1). In 
both the plots 8000 bins were used for averaging in the x-axis. ) 

The cumulated jump distributions show this very clearly. The cumulated jump length 
distributions in the chaotic structure regime show power-law behaviour (see Fig. [T2](a)). 
whereas the distributions in the mixing regime conform to a cumulated Gaussian (see Fig. 
[T2](b)). The same function fgaussi}) given above is used to obtain the cumulated distribution 
F{1). 

VI. TRANSPORT IN THE UDV REGION 

It was reported recently that the embedding map system studied here experiences a 
severe form of nonhyperbolicity in the aerosol regime in the neighbourhood of the parameter 
values a = 0.8, 7 = 0.40 [381]. A crisis is seen in the system at these parameter values. 
The nonhyperbolicity in higher dimensional systems, in general, manifests itself as unstable 
dimension variability (UDV) with an accompanying breakdown of the shadowing theorem. 
The signatures of the presence of unstable dimension variability in the embedding map were 
seen in the fluctuation of the Lyapunov exponent around zero. The presence of UDV was 
confirmed from the distribution of the finite time Lyapunov exponents which showed that 
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FIG. 12: The log- log plots of the cumulative probability distributions of the normalized jump 
lengths in (a) a typical trapping region (a = 0.5, 7 = 0.5) showing a power law tail with slope - 
2.5 and (b) a typical normal diffusion region (a = 2.8, 7 = 0.1) showing the fit with a cumulative 
gaussian distribution. In both the figures, 8000 bins were used for averaging. 



the spectrum is equally distributed in both positive and negative exponents [38]. 

The UDV region is found at an interface between two dynamical regimes, for example in 
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FIG. 13: The plot of A 

max vs 7 showing sudden jump at crisis. The inset shows the fluctuation of 
Amax around zero which is a signature of presence of UDV (fixed a = 0.8). 
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the above case, it is in seen at the interface between the periodic and the chaotic structure 
regimes. The transport properties of the system are studied in this region. However, due 
to the breakdown of shadowing at long times, computations are carried out for short times. 
Since the region where the UDV is seen lies between the periodic and chaotic structure 
regime, the statistical properties such as the recurrence times statistics and the diffusion 
properties will have contributions from both the periodic and the chaotic structure regime. 

The cumulative recurrence time statistics for a typical point in the UDV regime is shown 
in Fig. [MT a). There is an exponential decay of recurrence times at intermediate time scales, 
as seen for the recurrence time distribution in the chaotic structure regime. The periodic 
regime has its signature at the short and long time scales with the periodic behaviour showing 
up as discrete steps in the plot. The plot of the variance of the particle cloud as a function 
of time for the same parameters is seen in Fig. [T^ b). This is clearly a subdiffusive case 
with the slope rj ~ 0.6, indicating the influence of the chaotic structures. 




FIG. 14: (a) The cumulative recurrence time distribution showing an exponential fall in the inter- 
mediate time range. The discrete plateaus like behaviour in the small and large time scales are due 
to the periodic nature, (b) The variance as a function of time in the UDV case shows subdifFusvie 
behaviour. The parameter values are a = 0.8, 7 = 0.4062. 
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VII. CONNECTION BETWEEN THE DYNAMICAL AND TRANSPORT PROP- 
ERTIES 



It is clear from the previous discussion that there is an intimate connection between the 
dynamical and statistical properties of the system. We summarize our inferences in this 
section. 

The phase diagram of the dynamical regimes of the embedding map, a paradigm for the 
dynamics of inertial particles in fluid flows, contains three distinct dynamical regimes, viz. 
periodic orbits, chaotic structure regimes and mixing behaviour. There is also a clear distinc- 
tion between the dynamical behaviour in the aerosol regime, and the dynamical behaviour 
in the bubble regime. Similar distinctions are seen in the recurrence time distributions, and 
the diffusive and drift phase diagrams of the system. 

The recurrence time statistics are useful in the regimes where no regular periodic be- 
haviour is seen, i.e., in the regions where the largest Lyapunov exponent is greater than 
zero. For the embedding map, these are the chaotic structure regime and the mixing regime. 
The cumulative recurrence time distribution for the chaotic structure regime shows an ex- 
ponential decay followed by a power law tail. The exponential behaviour at short times 
reflects the mixing background in the inhomogeneous phase space. The power law seen here 
is the hallmark of the sticky regions in the phase space. The diffusion studies reveal that the 
chaotic structure regime shows subdiffusive character on the average. The distribution of 
jump lengths in the transient regime has an envelope that conforms to a Levy distribution, 
indicating that while short jumps dominate the distribution, the number of long jumps is 
sufficiently large to contribute a power-law tail. 

The mixing regime, as expected, shows the exponential decay of the recurrence time 
distribution. The transport in this regime shows normal or brownian diffusion, where the 
variance grows linearly with time. Normal diffusive behaviour is predominant in the mixing 
regime on both the aerosol and bubble sides. The distribution of jump lengths in this regime 
has a Gaussian envelope. 

Two phase diagrams, viz., the diffusion and the drift phase diagram, were constructed 
in the a — 7 parameter space to compare the nature of diffusion and dynamics in all the 
dynamical regimes of the embedding map. It was found that three main types of diffusion, 
viz normal (brownian type), sub-diffusive and super-diffusive types are seen in the system. 
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As mentioned above, normal diffusion is the dominant behaviour in the mixing regime on 
both sides of the phase diagram. Super-diffusive regions are seen in the periodic regime of 
the aerosols, and in a part of the periodic regime of the bubbles. However, the periodic 
regime on the bubble side also contains a regime where the inertial particles are trapped 
and have zero drift. This regime has sub-diffusive transport with stationary states. Here 
the variance saturates to a constant value after an initial sublinear rise. This regime shows 



behaviour consistent with early studies of impurity behaviour [39|, where it was observed 
that bubbles were pushed towards the islands forming regions of preferential concentration. 

Sub-diffusive transport is also seen in the chaotic structure regime on the aerosol side 
of the phase diagram where trapped states are also seen. However, the trapped states are 
non- stationary over here. Here, the variance fluctuates about a constant value after the 
initial sublinear rise. On the other hand, normal diffusion is predominant in the chaotic 
structure regime on the bubble side of the phase diagram. Ballistic diffusion similar to 
that seen in many natural processes is also seen in the embedding map. The regions of 
ballistic diffusion in the embedding map lie on the boundaries of the periodic regimes. The 
embedding map also shows the existence of unstable dimension variability at the boundary 
between the periodic regime and the chaotic structure regime. Thus, the signature of both 
the regimes are seen in the recurrence time distribution. However, the transport behaviour 
is clearly sub-diffusive. 



VIII. CONCLUSIONS 



The transport properties of bailout embedding map, a paradigm to model the dynam- 
ics of inertial particles, are studied in this paper. The base fluid flow is assumed to be 
incompressible and modelled by an area preserving map, the standard map. The resulting 
embedding map is dissipative and captures the qualitative dynamics of both particles that 
are heavier than the fluid, the aerosols, and lighter than the fluid, the bubbles. The main 
dynamical regimes in the embedding map system are reviewed and the statistical character- 
izers of transport in each distinct dynamical regime are evaluated. These are the recurrence 
time statistics and the diffusion and drift properties of system. An intimate connection is 
observed between the dynamical and statistical properties of the system. 

The predominant dynamical regimes in the system are the periodic regime, the chaotic 



24 



structure regime and the mixing regime. These, together with the nature of the inertial par- 
ticles, viz. whether aerosol or bubble, influence the diffusion, drift and recurrence properties 
of the system. The mixing regimes of the system show normal diffusion and exponential 
decay of recurrence times indicating short term correlations in the system. Chaotic struc- 
ture regimes posses inhomogeneous sticky regimes in the phase space, and these contribute 
power-law tails to the recurrence time distributions and jump length distributions. Superdif- 
fusive behaviour is seen in the periodic regimes of both aerosols and bubbles. However, some 
trapping regimes are seen in both the aerosol and bubble cases, and these show sub-diffusive 
behaviour and low drifts, including zero drifts in the case of stationary states. 

Our results may have implications in the context of realistic applications. It was found 
that the average drift is close to zero in the trapping regions. This may have important 
consequences in practical contexts. For example, consider flows with reacting impurities. 
Chemical species that have low diffusion and low drift rates can get localized and concen- 
trated in certain regions. These can then react with other chemical species and cause their 
enhancement or depletion, with further consequences for the environment. 

Rapid transport, as in the case of superdiffusive case, can provide a mechanism by which 
impurities can easily access wider regions in the available space. This maybe desirable 
or undesirable, depending on the context. For example smoke emanating from a chimney 
containing particulate matter with fast diffusion and slow drift can spread the pollutants in 
the proximity of the source rather than carrying it off to farther places. On the other hand, 
this regime may be highly useful in the context of nutrients spreading through a fluid. 

The phase diagrams of the system indicate that a rich variety of dynamical and statistical 
regimes are available to inertial particles in fluids. Specific choices of regime may be suitable 
for a specific application. Hence the insights gained from these diagrams may be useful in 
varied application contexts. 
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